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Abstract 

We have studied the synchronization in finite iV-unit FitzHugh-Nagumo neuron 
ensembles subjected to additive and multiplicative noises, by using the augmented 
moment method (AMM) which is reformulated with the use of the Fokker-Planck 
equation. It has been shown that for diffusive couplings, the synchronization may 
be enhanced by multiplicative noises while additive noises are detrimental to the 
synchronization. In contrast, for sigmoid coupling, both additive and multiplicative 
noises deteriorate the synchronization. The synchronization depends not only on 
the type of noises but also on the kind of couplings. 
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1 INTRODUCTION 



Nonlinear stochastic equations subjected to additive and/or multiplicative noises have 
been widely adopted for a study on real systems in physics, biology, chemistry, economy 
and networks. Interesting phenomena caused by both the noises have been intensively 
investigated (for a recent review, see Ref. 1, related references therein). It has been 
realized that the property of multiplicative noises is different from that of additive noises 
in some respects as follows. (1) Multiplicative noises induce the phase transition, creating 
an ordered state, while additive noises are against the ordering [2] -[6]. (2) Although 
the stochastic resonance is not realized in linear systems with additive noises, it may 
be possible with multiplicative color noise (but not with multiplicative white noise) [7, 
8]. (3) Although the probability distribution in stochastic systems subjected to additive 
Gaussian noise follows the Gaussian, it is not the case for multiplicative Gaussian noises 
which generally yield non-Gaussian distribution [9]-[12]. (4) The scaling relation of the 
effective strength for additive noise given by f3(N) = (3(1) /y/N is not applicable to that for 
multiplicative noise: a(N) ^ a(l)/y/~N, where a(N) and f3(N) denote effective strengths 
of multiplicative and additive noises, respectively, in the TV-unit system [13]. 

In order to show the above item (4), the present author has adopted the augmented 
moment method (AMM) in a recent paper [13]. The AMM was originally developed 
by expanding variables around their mean values in order to obtain the second-order 
moments both for local and global variables in stochastic systems [14]. The AMM has 
been successfully applied to a study on dynamics of coupled stochastic systems described 
by Langevin, FitzHugh-Nagumo and Hodgkin-Huxley models subjected only to additive 
noises with global, local or small-world couplings (with and without transmission delays) 
[15]. In Ref. [13], we have reformulated the AMM with the use of the Fokker-Planck 
equation (FPE), in order to avoid the difficulty due to the Ito versus Stratonovich calculus 
inherent for multiplicative noise. It has been pointed our that a naive approximation of the 
scaling relation for multiplicative noise: ct(N) = ct(l)/ VN, as adopted by Munoz, Colaiori 
and Castellano in their recent paper [6] , leads to the result which violates the central-limit 
theorem and which is in disagreement with those of AMM and direct simulations. 

The purpose of the present paper is two folds: (1) to reformulate AMM forFitzHugh- 
Nagumo (FN) model 

subjected to both additive and multiplicative noises with the use of FPE [13], and (2) 
to discuss the respective roles of the two noises on the synchronization. Our calculations 



2 



have shown that multiplicative noises may enhance the synchronization while additive 
noises work to destroy it. This is similar to the property in item (1) discussed above. 

The paper is organized as follows. In Sec. II, we have applied the DMA to finite N- 
unit FN networks subjected to additive and multiplicative noises. Numerical calculations 
are presented in Sec. III. The final Sec. IV is devoted to conclusion and discussion. 

2 Noisy FN neuron ensembles 
2.1 Augmented moment method 

We have adopted TV-unit FN neurons subjected to additive and multiplicative noises. 
Dynamics of a neuron i in a given FN neuron ensemble is described by the nonlinear 
differential equations (DEs) given by 

^ = F(x t )-cy t + aG(x t Mt)+f3Ut) + lt\t) + I {e \t), (1) 

% = bxi-dyi + e, (i = ltoN) (2) 

at 

with 

4 c) (*) = |Ete-*0- (3) 

In Eq. (l)-(3), F(x) = kx{x - a)(l — x), k — 0.5, a = 0.1, b = 0.015, d = 0.003 and 
e = [14] [16]: Xi and yi denote the fast (voltage) variable and slow (recovery) variable, 
respectively: G(x) an arbitrary function of x: I^(t) an external input whose explicit form 
will be shown shortly [Eq. (41)]: J expresses strengths of diffusive couplings, Z — N — 1: 
a and (3 denote magnitudes of multiplicative and additive noises, respectively, and rji(t) 
and express zero-mean Gaussian white noises with correlations given by 

(Vi(t) Vj (t')) = 5 ZJ 6(t-t>), (4) 

Ut) fc(0> = M(*-0, (5) 

(Vi(t)tjW)) = 0. (6) 

The Fokker-Planck equation p({xi}, {yi},t) is expressed in the Stratonovich represen- 
tation by[13][17] 



-gP = -^^{[F(x k )-cy k + I k ]p}-^—[(bx k -dy k + e)p] 
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where I k = + &\ 

We are interested also in dynamics of global variables X(t) and Y(t) denned by 

*(*) = ^E*i(*), (8) 

i 

nt) = ^E»(o- (9) 

i 

The probability of P(X, Y, t) is expressed in terms of p({ X i}, {yi}, t) by 

P(X,Y,t) = J Jdxdyi^iPi{xi^t)5{X-^Y, x i) 5 ( Y -^llyi)- ( 10 ) 

i i 

Moments of local and global variables are expressed by 

i x iVi) = =j J dx i dy y p i (x i ,y i ,t)x'ty-, (11) 

(X k Y £ ) = J J dXdYP{X,Y,t) X k Y l . (12) 

By using Eqs. (1), (2), (7) and (11), we get equations of motions for means, variances 
and covariances of local variables by 

^ = (F( Xl ))-c(y t ) + ^(G'( Xi )G(x t )), (13) 
d(y, 



dt 
dt 



= b{xi) - d{ Vi ) + e, (14) 



(xiF(xj)) + (xjF(xi)) - c(( Xi yj) + ( Xj yi)) 

J 
Z 



+ ^J2d X i X k) + (XjX k ) ~ ( x 1) - { x ])) 
k 

a 2 

+—[( Xi G\x i )G{x i )) + (xjG'ix^Gix,))} 



dt 

djxjyj) 

dt 



2 

+[a 2 (G(x i ) 2 )+(3 2 ]5 ij , (15) 

K&Vj) + ( x jVi)) - 2d (yiyj)> ( 16 ) 

(Vi F ( x i)) ~ c (yiVj) + H x i x j) - d(xi yj) 

2 

+| E(K Vj) - ( x i Vi)) + \(Vi G'( Xi )G(x t )), (17) 

where G'(x) = dG(x)/dx. 

Equations of motions for variances and covariances of global variables are obtainable 
from Eqs. (8), (9) and (12): 

~d~T = iv^ ( ^' (18) 

ijy^ = i EE f!<!!-M, (K , 7 = 1 , 2) (19) 

* 3 
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where we adopt a convention: Vu — Xi, Vn — Vi, V\ — X and V2 = Y. Equations (13) 
and (14) are used for N — 1 FN neuron (a = 0) and for N = 00 FN neuron ensembles 
(a = 0) in the mean-field approximation [18]. Equations (13)-(17) are employed in the 
moment method for a single FN neuron subjected to additive noises [16]. We will show 
that Eqs. (18) and (19) play important roles in discussing finite iV-unit FN ensembles. 
In the AMM [14] , we define the eight quantities given by 

= <K) = ^£<^>, (20) 

i 

7k,a = jjY,(( v Ki- Vk)(v m - (21) 

i 

p K ,x = ((V K -p K )(V x -p x )), («, 7 =1,2) (22) 

with 7i ; 2 = 72,1 and pi j2 = p 2 ,i. It is noted that 7 Kj a expresses the averaged fluctuations in 
local variables while p K>x denotes fluctuations in global variables. Expanding Eqs. (13)- 
(19) around means of p K as v Ki = \i K + 5v Ki , we get equations of motions for the eight 
quantities: 

^ - /o + /27i,i-^2 + ^ + /W (23) 



dt 
dt 



b/ii - dp 2 + e, (24) 



^ = 2(a 7lil -c7i, 2 ) + ^(pi,i-7i,i) + 2a 2 7i,i + « 2 /i? + /5 2 , (25) 
dt 



2(671,2-^72,2), (26) 



2(ap M - cp 1)2 ) + 2a 2 pi,i + + ^, (28) 
2(6pi, 2 -dp2 i2 ), (29) 



^7i 2 , / ,x JiV7 . a 2 7i 2 ,^„ N 

"g^ = 671,1 + («-rf)7l,2-C 7 2,2 + ^-(Pl,2- 71,2) + —^, (27) 

dpi 1 „. . „ 9 a 2 u? /9 2 

= 2(ap lil -c Pl , 2 ) + 2a 2 p lil + ^i + ^, 

dp 2 , 2 

U 1 f ^ 1 « 2 Pl,2 /Qn x 

-g- = + (a - a)Pi,2 - cp 2i2 + (30) 

where a = f\ +3/37^1, = and G(x) = x is adopted, relevant expressions 

for a general G(x) being given in the appendix. The original 2iV-dimensional stochastic 
equations given by Eqs. (1) and (2) are transformed to eight-dimensional deterministic 
equations. Equations (23)-(30) with a = (additive noises only) reduce to those obtained 
previously [14]. 
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2.2 N dependence of effective noise strength 

Comparing the (3 2 term in d^i^/dt of Eq. (25) to that in dp ltl /dt of Eq. (28), we note 
that the effective strength of additive noise is scaled by 

(31) 



As for multiplicative noise, however, the situation is not so simple. A comparison between 
the a 2 terms in Eq. (27) and (30) yield the two kinds of scalings: 

ol 

a — > — =, for iii term, (32) 



a — > a, for 7 1;1 and p^i terms, (33) 

The relations given by Eqs. (31)-(33) hold also for d^y 1:2 /dt and dpi^/dt given by Eqs. 
(27) and (30). Thus the scaling behavior of the effective strength of multiplicative noises 
is quite different from that of additive noises, as previously pointed out for Langevin 
model [13]. 

Nevertheless, we note that in the limit of J = 0, AMM equations lead to 

P«,a = («,A = 1,2) (34) 

which is nothing but the central-limit theorem describing the relation between fluctuations 
in local and average variables. 

2.3 Synchronization ratio 

In order to quantitatively discuss the synchronization, we first consider the quantity given 
by [14] 

= 4 £<M*) - ^(01 2 ) = 2[ 7 i,i(0 - PiM- (35) 

ij 

When all neurons are in the completely synchronous state, we get Xi(t) = X(t) for all 
i, and then R(t) = in Eq. (35). On the contrary, in the asynchronous state, we get 
R(t) = 2(1 — l/N)ji t i = Ro(t) from Eq. (34). We have defined the synchronization ratio 
given by [14] 

which is and 1 for completely asynchronous (R = Rq) and synchronous states (R — 0), 
respectively. We have studied the synchronization ratios at tf and t m as given by 

S f = S(t f ), (37) 
S m = S(t m ), (38) 
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with 



t f = {t | X{t) = 9,dX(t)/dt > 0}, (39) 
t m = {t | dS(t)/dt = 0}, (40) 

tj denoting the firing time at which the global variable X(i) crosses the threshold 9 from 
below and t rn the time when S(t) has the maximum value. Sf and S m depend on model 
parameters such as the noise intensities (a and (3), the coupling strength (J) and the size 
of cluster (N). 

3 CALCULATED RESULTS 

We have made numerical calculations, applying an external input given by 

&\t) = AQ(t-t rn )Q(t m + t w -t), (41) 

where A = 0.1, t in = 40 and t w = 10 [14]. AMM equations given by Eqs. (23)-(30) 
have been solved by using the fourth-order Runge-Kutta method with a time step of 
0.01. Direct simulations (DS) for the iV-unit FN model given by Eqs. (l)-(3) have been 
performed by using the Heun method with a time step of 0.003. Results of DS are averaged 
over 100 trials. All quantities are dimensionless. 

Figures l(a)-(d) show time courses of /ii(t), 7i,i(£), and S(t), respectively, 

calculated by AMM (solid curves) and DS (dashed curves) with a = 0.01, (3 = 0.001, 
J = 1.0 and iV = 100. When the external input J (e) (t) is applied at t = 40 to the 
quiescent states which have been randomized by applied small additive noises of (3 = 0.001, 
FN neurons fire, and 7i,i(t), Pi,i(t) and S(t) develop. Results calculated by AMM are in 
good agreement with those of DS. 

Figure 1(d) shows that when an input signal is applied at t = 40, S(t) is suddenly 
decreased but has a peak at t ~ 60 where X(t) is in the refractory period. In order to see 
the behavior of S(t) in more detail, we show in Fig.2, its time course for the four cases: 
(1) a = 0.0, (2) a = 0.002, (3) a = 0.01 and (4) a = 0.05 with (3 = 0.001, J = 1.0 
and iV = 100. In the case (1), the system is subjected only to additive noise, for which 
S(t) plotted by dashed curve is increased by an applied input for 40 < t < 50. It shows 
Sf = 0.30 at tf = 44.5 and S rn = 0.44 at t = 60.35, and approaches the equilibrium value 
of S = 0.159 at t > 100. In the case (2) with a = 0.002, S(t) shown by dotted curve yields 
Sf = 0.205 at t f = 44.5 and S m = 0.526 at t = 60.37. In the case (3) with a = 0.01, S(t) 
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shown by solid curve leads to Sf = 0.05 at tf — 44.5 and S m = 0.838 at t — 60.55. In the 
case (4) with stronger multiplicative noise of a = 0.05, S(t) plotted by chain curve yields 
Sf = 0.03 at tf — 44.5 and S m = 0.910 at t — 60.6. We note that with more increasing 
a, Sf is much decreased while S m is much increased. 

This trend is more clearly seen in Fig. 3(a), where Sf and S m are plotted as a function 
of a for (3 = 0.001, 0.01 and 0.02. In the case of f3 — 0.001, Sf is rapidly decreased and 
S m is rapidly increased with increasing a. For stronger additive noises of (3 — 0.01 and 
(3 = 0.02, Sf (S m ) is gradually decreased (increased) with increasing a. These results 
show that multiplicative noises enhance S m but deteriorate Sf. 

Figure 3(b) shows Sf and S m as a function of (3 for a = 0.0, a = 0.01 and a = 0.05. 
In the case of a — 0.0, Sf and S m are almost independent of (3 though they are gradually 
decreased for larger (3. In the cases of a = 0.01 and a = 0.05, Sf and S m are increased 
and decreased, respectively, with increasing (3. Although these results give an impression 
that additive noises enhance Sf and deteriorate S m , it is not true. Rather additive noises 
work to recover Sf and S m to the values of Sf = 0.30 and S m = 0.44 for the absence of 
multiplicative noise (a — 0). 

We have studied the effects of the noises on Sf and S m : the former expresses the 
synchronization ratio at the firing time at X(t) = 6 and the latter denotes its maximum 
value when X(t) is in the refractory period at t ~ 60. We may note that if multiplicative 
noises exits, the synchronization ratio S(t) is once decreased when an input applied, 
and it soon rebounds, showing the enhanced value. This trend is more significant for a 
considerable multiplicative noises. In this sense, the synchronization may be enhanced by 
multiplicative noises. 

By using our AMM, it is possible to study the dependence of the synchronization on 
the size of ensembles (N). Figure 4 shows the N dependences of Sf and S m in the two 
cases for (1) a = 0.01 and (3 = 0.001 (multiplicative noise dominant) and (2) a = 0.0 
and (3 = 0.01 (additive noise only) with J = 1.0. Both Sf and S m are increased with 
decreasing N, which shows that the synchronization becomes better for smaller systems. 

4 CONCLUSION AND DISCUSSION 

Although we have adopted the diffusive coupling given by Eq. (3), the sigmoid coupling 
given by 

4 C) W = f (42) 
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has been widely employed for discussing networks, where K expresses a coupling strength 
and H(x) is an arbitrary function of x [Eq. (51)]. Diffusive and sigmoid couplings model 
electrical and chemical synapses, respectively, in neuronal systems. It is worthwhile to 
study also the case of sigmoid coupling although such FN ensembles subjected only to 
additive noises were studied with the use of AMM [14]. 

A straightforward calculation using the AMM discussed in Sec. II leads to equations 
of motion given by 

' //f| = /o + Mi,i ~c/i 2 + i^ + />27i,i) + ^ + / (e) , (43) 



dt 
dp 2 



= bfj.1 - dp 2 + e, (44) 



2(a 7 i,i-c 7 i, 2 ) + ^i(pi,i-^) 

2_, i „.2,.2 i o2 



dj. 



2,2 



dt 



+2a 2 7l ,i + <* \A + P , (45) 
2(671,2-^72,2), (46) 



= 67i,i + (a-d)7i,2-c72,2 + ^i(/»i,2-^) 

(47) 

^ = 2{ap 1 , 1 -cp 1 , 2 ) + 2Kh lPl , 1 + 2a 2 p 1 , 1 + y^ + ^, (48) 

dp2,2 

dt 



= 2(b Pl , 2 -dp2, 2 ), (49) 



= + (a - d)/0i,2 - cp 2 ,2 + Kh,!p li2 + - , (50) 

where ^ = (1/! Comparing Eqs. (43)-(50) to Eqs. (23)-(30), we note the 

following points in coupling contributions between the two types of couplings: (i) the 
contributions in d'ji^/dt and d'ji^/dt terms of diffusive coupling are proportional to (pi,i — 
7i ; i) while those for sigmoid couplings are proportional to (pi ; i — 71,1/iV), (ii) dpi/dt, 
dpi,i/dt and dp^ 2 jdt terms in diffusive coupling have no contributions from the couplings 
in contrast to those in sigmoid coupling, and (iii) there are no differences in dp 2 /dt, 
d r ) 2 , 2 /dt and p 2)2 jdt terms. The item (i) mainly yields the difference between the effects 
of multiplicative noises on the synchronization for the diffusive and sigmoid couplings, as 
will be discussed shortly. 

We have performed numerical calculations by using Eqs. (43)- (50) with 

E[X) = {l + exp[(a;-0)H}' (51) 
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9 and w denoting the threshold and width, respectively. Time courses of /xi(t), 7i,i(t) 
and for 9 = 0.5, w = 0.1, K = 0.1 and iV = 10 are similar to those shown in Figs. 

1 and 6 in Ref. [14]. Figure 5 shows time courses of S(t) for three a values of a = 0.0, 
0.01 and 0.05 with (3 = 0.001. We note that S(t) has two peaks: one after an input is 
applied and the other when X(t) is in the refractory period. In the case of a — 0.0, we 
get S f = 0.108 at t = 44.16 and S m = 0.342 at t = 62.92. In the case of a = 0.01, we 
get S f = 0.073 at t = 44.16 and S m = 0.287 at t = 64.35. In the case of a = 0.05, we 
get S f = 0.053 at t = 44.15 and S m = 0.284 at t = 64.32. Both S f and S m are decreased 
with increasing a. This behavior is more clearly shown in Fig. 6 where S f and S m are 
plotted as a function of a. 

A comparison between Figs. 3(a) and 6 shows that with increasing a, S m for diffusive 
couplings is increased while that for sigmoid couplings is decreased. This difference mainly 
arises from the item (i) for the dji^/dt term discussed above. Namely, 7^1 for diffusive 
coupling is reduced by a negative contribution proportional to (pi,i —71,1) (i.e. pi t i < 7^1) 
which yields an enhancement in S(t) of Eq. (36). On the contrary, 7^1 for sigmoid cou- 
pling is slightly increased by a positive contribution proportional to (pi,i — 71,1/iV) (i.e. 
Pi,i > 7i,i/A0 which reduces S(t) of Eq. (36). In contrast, effects of multiplicative noise 
are not effective for Sf because Xj is not large at ij„ < t < tf. Then with increasing a, 
Sf is decreased for both the couplings. The main difference between the two couplings is 
the presence of the feedback (second) term of Eq. (3). Indeed, if we adopt the diffusive 
coupling without this term, which is equivalent to the sigmoid coupling with H(x) = x 
in Eq. (42), the synchronization is decreased with increasing a (result not shown). This 
situation is similar to that of the synchronization in small-world networks [19]. It was 
shown in Ref. [15e] that when a small-world network is made by introducing randomness 
to a regular network, the synchronization in the small-world network with diffusive cou- 
plings may be increased while that with sigmoid couplings is decreased. Our calculation 
implies that synchronization depends not only on the type of noises but also on the kind of 
couplings. This may also suggest that an ordered state in the multiplicative noise-induced 
phase transition reported in Refs. [2]- [6], might partly owe the diffusive couplings em- 
ployed in these studies: multiplicative noises could not yield an ordered state with sigmoid 
couplings . 

In summary, we have studied the synchronization in FitzHugh-Nagumo neuronal en- 
sembles subjected to additive and multiplicative noises, by reformulating AMM with the 
use of FPE [13, 14]. The property of the two noises in FN neuron ensembles is summarized 
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as follows. 

(a) the scaling relation: a(N) = a(l)/\/~N is not hold for multiplicative noises although 
the relation: @(N) = (3(1) /y/N is valid for additive noises, 

(b) multiplicative noises may enhance the synchronization (S m ) for diffusive couplings 
though both the two noises are generally detrimental to it, and 

(c) for both the additive and multiplicative noises, the synchronization is more increased 
in smaller N systems. 

The item (a) supplements the result for Langevin model [13]. The item (b) is similar to 
the item (1) of an ordered state created by multiplicative noises [2] -[6] mentioned in the 
introduction. 

A disadvantage of our AMM is that its applicability is limited to weak-noise cases. For 
multiplicative Gaussian noises, the probability distribution become non-Gaussian yielding 
divergent second and higher moments for a large a, to which the AMM cannon be ap- 
plied. On the contrary, an advantage of the AMM is that we can easily discuss dynamical 
property of the finite TV-unit stochastic systems. We have solved the eight-dimensional 
ordinary differential equations for FitzHugh-Nagumo neuronal ensembles. In contrast, 
within direct simulation and the FPE approach, we have to solve the 2iV-dimensional 
stochastic equations and the (2N + l)-dimensional partial differential equations, respec- 
tively, which are much laborious than AMM. Our AMM may be applied to a wide class 
of coupled stochastic models subjected to additive and/or multiplicative noises. 
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Appendix: Equations of motions for a general G(x) 

Although Eqs. (23)-(30) express equations of motion for G(x) = x, we present the 
result for a general form of G(x): 

= fo + f2li,i-cfi 2 + y[#o#i + 3(#i# 2 + 0002)71,1] +I {e) , (52) 
= bjii - dji 2 + e, (53) 



dt 

= 2(071,1 - C7i >2 ) + —=-(pi,i - 7i,i) + 2a {g 1 + 2^2)71,1 



dt 
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+ a 2 9 2 o+P 2 , (54) 

^ = 2(671,2-^2,2), (55) 

^ = 671,1 + (a - d)7i,2 - C72.2 + ^-(Pi,2 - 71,2) + y(# 2 + 2^2)71,1, (56) 

^ = 2(ap 1 , 1 -c Pl , 2 ) + 2a 2 (( 7l 2 + 2^2)pi,i + ^ + ^, (57) 

^ = 2{bp 1>2 -dp 2 , 2 ), (58) 

^ = 6pi,i + (a-rf)pi,2-cp 2 ,2 + ^-(^ 2 + 2^2)pi,2, (59) 



where g e = (l/£ !)G w (//i). For G(:r) = x, we get g = pi, #1 = 1, and g 2 = g 3 = 0, with 
which Eqs. (52)-(59) reduce to Eqs. (23)-(30). 
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Figure 1: (color online). Time courses of (a) /ii(t), (b) 7i i i(t), (c) Pi,i(t) and (d) S(t) 
for a = 0.01, = 0.001 with diffusive coupling (DC) of J = 1.0 and N — 100, solid 
and dashed curves denoting results of AMM and direct simulations, respectively. At the 
bottom of (a), an input signal is plotted. Vertical scales of (b) and (c) are multiplied by 
factors of 10~ 4 and 10~ 5 , respectively. 



Figure 2: (color online). Time courses of the synchronization ratio S(t) for a = 0.0 
(dashed curve), a = 0.002 (dotted curve), a = 0.01 (solid curve) and a = 0.05 (chain 
curve) with diffusive couplings of J = 1.0, (5 — 0.001 and N = 100. 



Figure 3: (color online), (a) a dependences of Sf and S m for (3 = 0.001 (circles), (3 = 0.01 
(squares) and (5 = 0.02 (triangles), and (b) f3 dependences of Sf and S m for a = 0.0 
(circles), a = 0.01 (squares) and a = 0.05 (triangles) with diffusive coupling of J = 1.0 
and iV = 100: results of S m and S a of AMM are expressed by solid and chain curves, 
respectively, and those of DS by filled and open marks, respectively. 



Figure 4: (color online). iV dependences of Sf (solid curve) and S m (chain curve) for 
two sets of parameters: (1) a = 0.01 and j3 = 0.001 and (2) a = 0.0 and (3 = 0.01 with 
diffusive coupling of J = 1.0 and iV = 100, calculated by AMM and DS (circles and 
squares) . 



Figure 5: (color online). Time courses of the synchronization ratio S(t) for a = 0.0 
(dashed curve), a = 0.01 (solid curve) and a = 0.05 (chain curve) with sigmoid couplings 
(SC) ofK = 0.1, = 0.001 and N = 10. 



Figure 6: (color online), a dependences of Sf and S m for P = 0.001 (circles) and P = 0.01 
(squares) with sigmoid coupling of K = 0.1 and iV = 10: results of S m and S a of AMM 
are expressed by solid and chain curves, respectively, and those of DS by filled and open 
marks, respectively. 
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